function target_values=calculate_target_values(individual)
global problem
    switch problem
        case 'SCH'
            value1=individual(1)*individual(1);
            value2=(individual(1)-2)*(individual(1)-2);
            target_values=[value1,value2];
        case 'FON'
            sum1=0;
            sum2=0;
            for i=1:3
                sum1=sum1+(individual(i)-1/sqrt(3))^2;
                sum2=sum2+(individual(i)+1/sqrt(3))^2;
            end
            value1=1-exp(-sum1);
            value2=1-exp(-sum2);
            target_values=[value1,value2];
        case 'POL'
            A1=0.5*sin(1)-2*cos(1)+sin(2)-1.5*cos(2);
            A2=1.5*sin(1)-cos(1)+2*sin(2)-0.5*cos(2);
            B1=0.5*sin(individual(1))-2*cos(individual(1))+sin(individual(2))-1.5*cos(individual(2));
            B2=1.5*sin(individual(1))-cos(individual(1))+2*sin(individual(2))-0.5*cos(individual(2));
            value1=1+(A1-B1)^2+(A2-B2)^2;
            value2=(individual(1)+3)^2+(individual(2)+1)^2;
            target_values=[value1,value2];
            
        case 'KUR'
            value1=0;
            value2=0;
            for i=1:2
                value1=value1+(-10*exp(-0.2*sqrt(individual(i)^2+individual(i+1)^2)));
              
            end
            for i=1:3
                value2=value2+abs(individual(i))^0.8+5*sin(individual(i)^3);
            end
            
            target_values=[value1,value2];
            
        case 'ZDT1'
            n=30;
            value1=individual(1);
            
            
            g=1+9*(sum(individual)-individual(1))/(n-1);
            value2=g*(1-sqrt(individual(1)/g));           
            target_values=[value1,value2];
        case 'ZDT2'
            n=30;
            y=individual(2:n);
            gx=1+9*sum(y)/(n-1);
            value1=individual(1);
            value2=gx*(1-(individual(1)/gx)^2);
            target_values=[value1,value2];
            
        case 'ZDT3'
            %ZDT3
            n=30;
            y=individual(2:n);
            gx=1+9*sum(y)/(n-1);
            value1=x(1);
            value2=gx*(1-sqrt(individual(1)/gx)-individual(1)/gx*sin(10*pi*individual(1)));
             target_values=[value1,value2];
            
        case 'ZDT4'
            %ZDT4
            n=10;
            y=individual(2:n);
            gx=1+10*(n-1)+sum(y.*y-10*cos(4*pi*y));
            value1 = individual(1);
            value2 = gx*(1-sqrt(individual(1)/gx));
            target_values=[value1,value2];
    end
    

end